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Electron density profiles in the afterglow of a stabilized linear pinch have been measured using 
the refractive properties of the plasma column. The deflection angle of a microwave beam refracted 
by the discharge is observed as a function of beam position. The resulting data lead to an integral 
equation for the electron density distribution, which is reducible to the Abel integral equation. The 
electron density is derived from the experimental data by inverting this last equation. 

The method has been checked using direct transmission and interferometric measurements. 



1. Introduction 

Various branches of plasma physics, as well as 
the diverse applications of this science, call for a 
detailed knowledge of the electron density distribu- 
tions occurring in a number of experimental situations. 

Electromagnetic waves in a plasma are affected 
by the electrons present in it, so that microwave and 
optical diagnostic techniqes lead to data about elec- 
tron densities. Present day sources in this part of 
the electromagnetic spectrum are characterized by 
a relatively narrow band, which prohibits the use of 
well-known ionospheric techniques employing large 
frequency excursions in the 10 12 — 10 16 el/cm 3 density 
range [Budden, 1961]. 

The problem of electron density profile measure- 
ment at a fixed frequency has been attacked in the 
past both theoretically and experimentally. The 
method described by Kharadly [1963] was developed 
for a plane stratified medium, and consists in measur- 
ing the distance between the incident and emergent 
beams as a function of the angle between the incident 
beam and the normal to the strata (fig. la). The 
method can also be applied to a cylindrical discharge, 
the ray path in figure la being contained in a meridian 
plane. The method can be shown to lead to unique 
electron density profiles, provided the profile is a 
monotonic function of position. 

The method proposed by Shmoys [1961] makes 
use of the refractive properties of the plasma, the 
deflection angle 6 being recorded as a function of the 
offset distance b between the incident ray and the 
midplane of the discharge. The electron density 
distribution can be deduced uniquely from the 6(b) 
dependence (fig. lb), again provided the profile is 
a monotonic function of position. A common char- 
acteristic of both methods is that the density profile 
is deduced from a set of measured geometrical quanti- 
ties at a single microwave frequency. In both methods 
it is necessary to make assumptions about the nature 



of the electron distribution. Working in the meridian 
plane, one must assume that the density does not 
change along the strata and this requires the use of a 
sufficiently long tube. Axial symmetry is assumed 
in the transverse plane method, and measures have 
to be taken to ascertain the correctness of this assump- 
tion. 

A similar approach was adopted by Wort [1963] 
who limited himself to the case of a parabolic electron 
density and applied this theory to the diagnostics of a 
statistically distributed bundle of plasma columns in 
a turbulent plasma. 

This paper deals with the application of the trans- 
verse plane method to the measurements of electron 
density profiles in the afterglow of a linear pinch 
discharge. Section 2 deals with the theory of the 
experiment i.e., the calculation of density profiles 
from refraction data. It is essentially the same 
theory as that by Keller et al. [1956], except for the 
fact that the number of integrations in the final result 
has been reduced from two to one. The experiment 
itself is described in section 3. 
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FIGURE 1. (a) Plane stratified medium. 
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2. Theory 

The ray trajectory is derived from Snell's law which 
has the form 



I(r) r sin <p = const. 



(1) 



in cylindrical coordinates; I(r) is the refractive index, 
assumed to be a function of r only, <p is the angle 
between the ray trajectory and the radial direction. 

Now 



r dO 



Ctgip 



(2) 



which, combined with (1), gives the differential equa- 
tion of the trajectory. We assume a collisionless 
plasma, so that the index of refraction can be expressed 
through the electron density only 



P(r) = l 



n 



(3) 



where n c is the critical electron density corresponding 
to the frequency of the microwave signal 

n c = meo(A) 2 le 2 . 

Combining (1), (2), and (3) and integrating, we get 
the expression for the deflection angle 6(b) (fig. lb) 



0(b)- 



ri <& 



(4) 



which is derived by Shmoys [1961] . a is the radius of 
of the plasma column, i.e., that radius at which the 
electron density drops to zero; b is the offset distance 

defined in figure lb; as <^ = sin _1 — is the initial angle 

r 
between the ray and radius vector (at r = a); x = - 

a 

is the normalized radius, xo is the normalized radius 

at the point of closest approach; i.e., the smallest 

real root of 



1 



n(x) 
n c 



■1 = 0. 



(5) 



Equation (4) is an integral equation involving the elec- 



If the observed quantities are a&fao) and <p , the 

two experiments become identical for large values of a. 

Equations (6) and (7) enable us to bring (4) in the form 



*(6) 



dx 



2 sin <p jx 1 2 (, n\ 2 

x -dx 2 1 1 1 — 1 + cos z <p 
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ntroducing the auxiliary variables 




Y- * (6) . 

2 sin <po 


(9) 


U=cos 2 y?o> 


(10) 


r-i-*(i-3 


(11) 



we bring (4) to the standard Abel form 

1 dx 
~~;dV 



Y(U) 



= _ f" xdV 

Jo vzrn 



:dV. 



(12) 



The lower limit in the last equation is zero as 

n(l)\ 



F(l) = l-1 2 1 



= 0, 



since we assume n(l) = 0. If n(l) is finite, (4) must 
be modified (invoking Snell's law) to take into account 
the finite refraction at the plasma — free space 
boundary. 

Equation (12) is of the Abel type and has the inver- 
sion [Hamel, 1949] 

1 dx = _l__d_ [ w Y(U) 

xdW ndWJo VW F ^U (13) 

W is a parameter to be substituted by V after integra- 
tion and differentiation. The last equation can be 
integrated between the limits and W: 



log*(JF)-logs(F=0) = " f l 

7T Jo 

In a penetrable plasma 



Y(U) 

Vw-u 



dU. 



and 



1 



n c 



< 1 



can be 1 onlv if x- 



F=l-x 2 (l 
1 and n = 0. 



(14) 



n 
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In an 



tron density distribution. In the case of the plane impenetrable plasma the same applies to the sub 



st atified medium, one is also led to an integral equa- 
tion for the density profile which can be reduced to 
an equation of the Abel type. To increase the simi- 
larity between the two (mathematically different) 
cases, we introduce the variables 



critical region where — < 1. For the overcritical 
region 1 <0 and V is always positive. So the 

TLq 

second term — log x(V=0) at the left-hand side of 
(14) is always zero, which gives 



*(&) = 7r-2<p<>-0(&)> 

b 

sin <£o — — * 



(6) 
(7) 
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x{W) = e\p 



_1 f» 

"■Jo 



Y(U) 



Vw-u 



dU 



(15) 



Note that only one integration has to be performed 
in (15). It gives the answer to the problem of finding 
the density distribution from refraction data given in 
the form of a Y(U) function; x(V) can be computed 
from (15) taking arbitrary values of V and computing 

the corresponding values of — from (11). 

n c 

Let us check (15) in the simplest case of no electron 

density. It is evident that 6(b) = and consequently 

^ = 7r— 2<£> which gives 



Y(U)-- 



V 



77 — 2y>p _77" 

sin ^° 2Vl-cosVo~ 



•2 arc cos vU 



2V1-U 



Y(U) = 



arc s 



mW 



Vi-u 



(16) 



As is measured as a function of b we obtain another 

useful expression by introducing 5 = sin <fo = ~ instead 

a 

of U=cos 2 <£o as the variable. To simplify the result 
we put s = Vl — i/, w = Vl — W and get 



(21) 
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He 


x 2 — v 2 

X 2 





(22) 



It should be observed that the integration in (21) 
is elementary for a polynominal 9{s). The actual 
procedure of interpreting the experimental data in- 
volved a polynomial of second degree. 



The x{W) dependence is obtained from (15): 
arc sin \U 



x(W) = exp 



1 f w 
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V(l-U)(W-U) 



dU 



Now, it can be proved that 
" w arc sin VZ/ 



Jo 



(17) 



-</{/=-- log (l-r) 

^/(\-U)(W-U) 2 



(18) 



giving, after substitution of W by V 



x{V) = yJ\-V ■ 



3. Experiments 

The basic idea of the experiment is shown in figure 
lb. Figure 2 is a schematic diagram of the microwave 
setup. A narrow microwave beam from the trans- 
mitting antenna is refracted by the plasma column and 
picked up by the receiving antenna. The trans- 
mitting horn is mounted on a microwave bench allowing 
lateral displacement of the incident beam. The 
receiving horn slides on another microwave bench and 
is always displaced for the same distance b as the 
transmitting horn. This second bench is fastened to a 
turntable fitted with an angular dial permitting the 
reading of the angle 0. 



The electron density is found from (11): 

v=\-a-v)(i-^\ 

i.e., n = 0. 

The general equation (15) can be simplified using the 
integral (18). From (6) and (15) we get 



x{W) = exp \ 



1 



-2y?o 



^Jo V(i-u)(w-u) 

J_ f w 6(b) 

277J0 



:dU 



y/a- U)(w-U) 



dU\ 



(19) 



The first integral is the result (18), and we get 

6(b) 



Mm -VT=We* P {±l w w 



■U){W-U) 



dU 



(20) 




Figure 2. Microwave setup. 

A — attenuator. 

I — ferrite isolator. 

K - klystron. 

D — directional coupler. 

W — wavemeter. 

C — crystal. 

S — short circuit. 

TA — transmitting antenna. 

RA — receiving antenna. 

CRO — cathode ray oscilloscope. 
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Figure 3. Linear pinch equipment. 



rectifier unit. 
Helmholtz coils, 
shunt, 
capacitance divider. 




Figure 4. Setup for checking axial symmetry. 




Figure 5. Check of axial symmetry. 



The signal of the right-hand side horn is negative. The time scale is 0.5 
msec/cm, the voltage scale is 0.05 V/cm. 



The plasma was produced by discharging a con- 
denser bank (180 jxF, 6 kV max.) through argon at 
75 mtorr pressure, and hydrogen at 107 mtorr. A 
spherical vessel 70 cm in diameter fitted with pointed 
electrodes was used. The electrodes were fed through 
a large inductive loop resulting in a damped oscillatory 
current waveform, with a period of approximately 
200 /usee and peak currents of 10,000 A. An A.E.I. 
ignitron, type BK24A, was used as the switching 
element. 

It was found essential to use a stabilizing axial 
magnetic field during the pinch process. This field 
(fig. 3) was produced by passing 10 kA currents through 
a pair of single-turn Helmholtz coils, and had a maxi- 
mum value of 200 G. The formative time lag of the 
main discharge is of the order of 20 /usee during 
which period the axial magnetic field is allowed to 
increase, being trapped at the start of the main dis- 
charge. Without this precaution the plasma in the 
afterglow was found to be nonreproducible and axially 
nonsymmetric. As axial magnetic fields are known 
to suppress instabilities in linear pinches, the magnetic 
field is considered to stabilize the initial conditions of 
the discharge, leading to an axially symmetric plasma 
in the afterglow. 

As the main assumption in the theory of the experi- 
ment is the axial symmetry of the plasma column, this 
was checked first. The arrangement is shown in 
figure 4. Two receiving horns at 45° from the direction 
of the transmitting horn pick up the refracted micro- 
wave signals. It was found that both left and right 
signals occurred at approximately the same time 
(fig. 5), which ascertains axial symmetry at that parti- 
cular instant of time. 

Another essential assumption is the reproduci- 
bility of the discharge. Values of 6 corresponding to 
different values of the offset distance b are measured 
in different firings, so that consistent results can be 
expected only if the electron density is sufficiently 
reproducible from one discharge to another. Repro- 
ducibility was constantly checked by observing the 
pulse shape and the peak of pulse time in the sym- 
metry experiment. Even with an axial magnetic 
field the first few discharges were dissimilar to all 
other subsequent firings. Discharges took place at 
regular time intervals (1 min). After at most 10 dis- 
charges, the pulses did not change in shape and time. 
If the apparatus was not operated for 15 min, this 
effect was not observed, but an interval of 1 hr was 
sufficient to cause a few initial pulses to be irregular. 
This effect was ascribed to the gases occluded in the 
electrodes. 

All experiments have been carried out at the fixed 
wavelength of 8.54 mm. Both transmitting and 
receiving aerials were phase-corrected horn-lens 
combinations, designed on the basis of geometrical 
optics. The lens diameter was 70 mm. A lens 
diameter of the order of 10X is essential, as smaller 
lenses yield unreasonably diverging beams. 

Measurements at a prescribed time in the afterglow 
can be performed by observing the angle as a func- 
tion of time with the offset distance b as a variable 
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parameter and then determining 6 as a function of b 
at a given time from the resulting curves. A different 
procedure was adopted here however. Both receiving 
and transmitting horns arc first set to a given offset 
distance b. Then the turntable is manipulated in the 
course of firing until the peak of the received microwave 
pulse falls at the required instant of time on a cathode 
ray oscilloscope. The sensitivity of the method 
depends upon the narrowness of the received pulses. 
The pulse shape changes appreciably with time in the 
afterglow and also with the bank voltage. Best results 
in argon have been obtained at 1600 ^isec after the 
start of the discharge with a bank voltage of 3.5 kV. 
Direct transmission experiments indicate that the 
discharge was critical at 1500 /xsec. 

The measured 0(b) curve is the solid line in figure 6. 
Note that there is no refraction for b/a > 0.35. The 
curve is the mean of 5 measurements, the maximum 
deviation in the results being 3°. To obtain the elec- 
tron density profile a parabola of the type 



0=0o 



1 



5 < 5 , S : 



b 

a 



6 = 0, s > s 

was fitted to the experimental curve by least square 
methods, giving 



0o = 44° 



5o = 0.347. 



The dotted line in figure 6 is the parabolic approxi- 
mation. Equations (21) and (22) have been used to 
compute the electron density profile, giving 



x(v) =v exp 



77 



v^ 



2 5o 






n 
n c 



The profile obtained from the last two equations is 
the solid line in figure 7. The fact that n = n c at r = 
is a result of the failure of geometrical optics near the 
axis of the plasma column. In ray theory any 0(b) 
curve with a finite value of 0(0) leads to n = n c at the 
axis. 

The anomaly may be understood by referring to 
figure 5 in Shmoys [1961]. This is a 6(b) curve for a 
parabolic distribution function, and it is fairly reason- 
able to assume that the electron density near the center 
of the plasma column is parabolic. We see from the 
curve that most of the electromagnetic energy in the 
central part of a finite beam would be scattered be- 
tween = and 6 = 6 max , the outer portion of the beam 
forming a caustics at 9 = 6 m . dX . The maximum energy 
flow is not to be expected near 6 = 0, although the 
plasma is penetrable. 



~~^>! 






















"SX 












^ 












\ 

\ 



2.5 

0.07 



5 
0.14 



7.5 
0.21 



10 
0.28 



12.5 b(cm) 
0.35 b/a 



FIGURE 6. Measured 0(b) dependence. 

Argon pressure 75 mtorr, bank voltage 3.5 kV. The curve is the mean of 5 
measurements the maximum deviation in the results being 3°. Note that 
there is no refraction for b/a > 0.35. The dotted line is a parabolic approxi- 
mation of the experimental data obtained by least square methods, 

0«0o [1 -(«/*>)*], 0e = 44 ft , «o -0.347, **=£. 
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Figure 7. Electron density profile in the afterglow of a linear pinch 
in argon. 

The solid line is obtained by inverting the integral equation for the electron 
density distribution, assuming a parabolic 8(b) dependence (fig. 6). The 
point labeled IM is a result of interferometric measurements. The dotted 
line is a parabola, representing the density near the column center. 



The electron density at x = was derived from 
interferometric measurements. The interferometer 
shown in figure 2 is set to measure the phase shift of 
the microwave signal reflected from a plane mirror 
at b = 0. Figure 8 is a logarithmic plot of phase 
shift versus time. We deduce from this figure that 
the electron density at the axis was 0.8 n c 100 /xsec 
after the plasma was critical. This is the point labeled 
IM in figure 7. The broken line in the same figure 
is a parabolic continuation starting at that point. 

For b/a > 0.35 the plasma is not sufficiently refrac- 
tive to cause any measurable deflection of the beam. 
To verify this result an independent check of the 
experiment was provided. This is essentially a direct 
transmission test. The transmitting and receiving 
horns are located opposite each other on a line 
b = const., and the time during which the microwave 
signal is cut off is recorded. 
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FIGURE 8. Phase shift versus time from interferometric measurements. 

The average electron density, proportional to phase shift, drops for 20 per- 
cent in a 100 /usee interval. 
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FIGURE 9. 0(b) measured in hydrogen. 

Pressure 107 mtorr, bank voltage 4 kV. 
the curve is the parabolic approximation. 
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electron density as < — — n c at t = 1600 /usee from the 

decay rate obtained from interferometric measure- 
ments. The fact that the microwave beam at the edge 
of the vessel is unperturbed during the discharge and 
the whole of the afterglow provides an experimental 
test of the assumption /?(1) = 0. 

A similar experiment was carried out in hydrogen at 
107 mtorr. Unlike argon, where the visible discharge 
was confined to a narrow column, the light from the 
hydrogen discharge came from a central area occupy- 
ing 1/3 of the vessel. Without the axial magnetic 
field the luminous discharge fills practically the whole 
of the vessel, but the reproducibility of the microwave 
cutoff time was not sufficient to permit measurements 
of profiles. With the axial magnetic field, best repro- 
ducibility and symmetry have been obtained at 
4 kV condenser bank voltage. The plasma is critical 
at 245 jjisec. The 0(b) was measured at 250 jitsec 
after the application of the voltage pulse to the gap. 
The curve is shown in figure 9. The points are the 
mean of three measurements with a maximum devia- 
tion of 5°. The solid line is the parabolic approxi- 
mation given by the equation # = 71 — 90 s 2 . The 
profile obtained from this approximation using (23) 
and (24) is given in figure 10. For r=0, n = n c again, 
but there is no cusp at this point. It is a property 
of the parabolic approximation used that it leads to 
spiky profiles for 0(0) < 60°. Irrespective of the 
character of the derivative of the profile curve near the 
axis the value of the electron density obtained in this 
region are not correct unless the beam width is smaller 
than the radius corresponding to the maximum in the 
theoretical 6(b) curve. According to interferometric 
measurements the electron density at the axis is 
0.94 n c . This is the point labeled IM in figure 10. 
The direct transmission test applied to the hydrogen 
discharge shows that the electron density at the wall is 
zero. 

A suitable quantitative criterion for the applicability 
of ray optics is necessary both for the argon and hydro- 
gen discharges. Experimentally one measures the 
maxima of the received pulses. In ray theory, these 
coincide with the central ray of the beam. Wave 
optics gives different positions for the maxima. 
The quantitative corrections to ray theory necessary 
from the standpoint of wave theory have been com- 
puted by A. L. Cullen [1965] for a plane stratified 
medium. It follows from his theory that the applica- 
bility of ray optics to the oblique incidence experiment 
in critical plasmas of the sizes encountered here is 
justified. 



FIGURE 10. Electron density profile in the hydrogen discharge. 

The point labeled IM is the result of interferometric measurements. 

This was done at 6 = 5, 10, 15, 20, 25, and 30 cm. 
The plasma does not perturb the microwave beam at 
25 and 30 cm. At 15 and 20 cm the plasma was found 
penetrable for periods sufficiently long to estimate the 



4. Conclusions 

A microwave diagnostic method giving electron 
density profiles was applied to the afterglow plasma 
of a stabilized linear pinch. The deflection of a 8mm 
microwave beam is measured as a function of beam 
position. A simplified theory of the experiment is 
presented. 
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The ray theory version of the method is applicable 
only to very large, spatially slowly varying plasmas. 
With a finite microwave beam width difficulties occur 
near the center of the column where the deflection of 
a ray varies very much with ray offset distance. In a 
time varying plasma the method yields electron den- 
sity profiles only in a short interval of time in which 
the plasma is sufficiently refractive, in the vicinity of 
the instant at which the center of the column is critical. 
The sensitivity of the method is affected by the shape 
of the received microwave pulses and the reproduci- 
bility of the discharge. The procedure does not re- 
quire a great deal of axial space. On the other hand, 
the plasma column must possess symmetry of revolu- 
tion, and this assumption must be checked. 
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